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Abstract 

An irreversible cycle model of the Stirling refrigeration-cycle, using a ferroelectric material as 
the working substance, is established. Several irreversibilities due to thermal resistances 
between the working substance and the heat reservoirs, regenerative losses in two regenerative 
processes, and heat-leak loss between the heat reservoirs are taken into account. The influence 
of these irreversible losses on the performance of the ferroelectric Stirling refrigeration-cycle is 
analyzed, based on a general expression of the polarization of ferroelectric materials and a linear 
heat-transfer law. The cooling rate is optimized for a given power input. Some fundamental 
optimal relations and general performance characteristic curves of the cycle are obtained. The 
maximum cooling rate and other relevant performance parameters are determined. Some special 
cases are discussed in detail. When the regenerative losses are neglected, the results obtained 
may be directly used to describe the optimal performance of a ferroelectric Carnot refrigera¬ 
tion-cycle. Moreover, it is expounded that the calculated results are very general and also sui¬ 
table for the ferroelectric Stirling and Carnot refrigeration cycles whose working substances 
obey the Curie-Weiss law and Curie law. © 2002 Elsevier Science Ltd. All rights reserved. 


1. Introduction 

Since the early 1970s, many researchers have been searching for more environment- 
friendly refrigerants and new advanced cycle modes [1^1], Ferroelectric refrigerators 
use no environmentally sensitive materials and consequently have received increasing 
attention [5-10]. The previous investigations were mainly concentrated on the 
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Nomenclature 

a 

= Ql/ As, defined in Eq. (20) 

b 

= (1 — tj)Cpa/As, defined in Eq. (20) 

b\ 

= 2u/As, defined in Eq. (21) 

C 

Curie constant 

Cp 

specific heat at constant electrical polarization 

CpA 

average value of Cp 

d 

= 0(1 + b)(b + e)T c /T h - be, defined in Eq. (29) 

E 

electric-field strength 

Ei 

electric-field strength in the state i( = a, b, c, d) 

e 

= bjkT c , defined in Eq. (27) 

F(P) 

function determined by Eq. (1) 

h 

= \Jk\Jki 

k 

= 1/(1 /V^T + l/V^h) 2 ’ defined in Eq. (27) 

kn 

Boltzmann constant 

k L 

coefficient of heat leak 

k\ 

thermal conductance between the working substance and the heat 
sink 

k 2 

thermal conductance between the working substance and the cooled 
space 

N 

number of electric dipoles per unit volume 

P 

electrical polarization 

Pi 

electric polarization in the state i( = a, b, c, d) 

P 

power input 

Pr 

power input at maximum cooling rate 

Ps 

power input at zero coefficient-of-performance 

P* 

dimensionless power input 

Qc 

net heat absorbed from the cooled space per cycle 

Qh 

net heat released to the heat sink per cycle 

Ql 

heat leak loss per cycle 

Q, 

regenerative losses per cycle 

Qx 

heat released to the heat sink in the high-temperature isothermal 
process 

02 

heat absorbed from the cooled space in the low-temperature iso¬ 
thermal process 

012 

heat transferred into the regenerator per cycle 

034 

heat transferred out of the regenerator per cycle 

1L 

rate of the heat leak 

r 

cooling rate 

^max 

maximum cooling rate 

r* 

dimensionless cooling rate 


entropy of ferroelectric materials per unit volume 

^0 

entropy of ferroelectric materials per unit volume when P = 0 
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Si 

entropy of ferroelectric materials per unit volume in the state i( = a, 
b, c, d) 

As\ 

entropy change of ferroelectric materials per unit volume in the high- 
temperature isothermal process 

As 2 

entropy change of ferroelectric materials per unit volume in the low- 
temperature isothermal process 

T 

absolute temperature 

T c 

temperature of the cooled space 

T h 

temperature of the heat sink 

Co 

Curie temperature 

Ci 

temperature of the working substance in the high-temperature iso¬ 
thermal process 

C 2 

temperature of the working substance in the low-temperature iso¬ 
thermal process 

T\ r 

value of T\ at maximum cooling rate 

T lr 

value of T 2 at maximum cooling rate 

t 

cycle period 

h, hi hi U 

times spent on two isothermal and two regenerative processes 
respectively 

u 

a proportional constant appearing in Eq. (18) 

V 

Lagrange multiplier 

y 

= Cj/72 

P 

a parameter which is independent of temperature 

s 

coefficient of performance 

Sr 

coefficient of performance at maximum cooling rate 

>1 

efficiency of the regenerator 

h 

moment of electric dipole 


ferroelectric Carnot cycle. Like the magnetic Stirling refrigeration cycle, it is expected 
that the ferroelectric Stirling refrigeration-cycle will play an important role in 
refrigeration in the near future. 

The Stirling refrigeration-cycle is one of the important cycle models of refrig¬ 
erators with regeneration. According to the theory of classical thermodynamics, the 
Stirling refrigeration-cycle may possess the condition of perfect regeneration 
through the use of a reversible regenerator. Its coefficient-of-performance is equal to 
that of a reversible Carnot refrigeration-cycle. However, when the irreversibility of 
the heat transfer is considered, the performance of the Stirling refrigeration cycle is 
different from that of a Carnot refrigeration-cycle [11-13]. In recent years, although 
the influence of irreversibilities on the performance of the Stirling refrigeration cycle 
using a gas or magnetic material as the working substance has been investigated 
extensively, the influence of irreversibilities on the performance of the ferroelectric 
Stirling refrigeration-cycle has been rarely studied. 
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In the present paper, an irreversible cycle model of the ferroelectric Stirling refrig¬ 
erator will be established and used to investigate the influence of irreversibilities due 
to finite-rate heat transfers between the working substance and the heat reservoirs, 
regenerative losses in two regenerative processes, and heat-leak loss between the heat 
reservoirs on the performance of the cycle. The results obtained here will be of real 
significance. 


2. Thermodynamic properties of ferroelectric materials 


According to the theory of the Boltzmann statistics, the ferroelectric polarization 
may be given by [14,15] 


P — Njx tanh 


iAE + PP) 
k B T 


( 1 ) 


where N is the number of the electric dipoles per unit volume, /x is the moment of an 
electric dipole, k B is the Boltzmann constant, and is a constant which is indepen¬ 
dent of temperature, E is the electric-field strength, and T is the absolute temperature. 

Using Eq. (1), one can prove [10] that the entropy of ferroelectric materials per 
unit volume and the specific heat at constant electrical polarization are, respectively, 
given by 

5 = So (T) - P(E+ PP)/T+ M-lncosh + f P) (2) 

k B T 

and 

C p = Cp(T), (3) 

where so(T) is the entropy of ferroelectric materials per unit volume when P = 0 and 
is a function only of temperature. Its value depends on the properties of the specific 
ferroelectric material. Eq. (3) shows clearly that the specific heat at constant 
electrical polarization is a function only of temperature while independent of the 
properties of the ferroelectric materials. 


3. An irreversible ferroelectric Stirling refrigeration-cycle 

When a ferroelectric material is used as the working substance, the Stirling 
refrigeration-cycle consists of two isothermal and two constant electrical polariza¬ 
tion processes, as shown in Fig. 1, where 7), and T c are the temperatures of the heat 
sink and the cooled space, T\ and 73 are the temperatures of the working substance 
in the high- and low-temperature isothermal processes, Q\ and Qi are the amounts 
of heat exchanged between the working substance and the heat reservoirs during two 
isothermal processes, Qu and Q 34 are the amounts of heat exchanged between the 
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working substance and the regenerator during two constant electrical polarization 
processes. 

When the cycle is reversible, the temperatures T\ and T 2 of the working substance 
in two isothermal processes are identical with the temperatures 7), and T c of the heat 
reservoirs, respectively. Using Eq. (2), one can calculate that the amounts of heat 
released from and absorbed by the working substance during the high- and low- 
temperature isothermal processes are, respectively, given by 

Q\ = T\(s a - s h ) = TiAsi (4) 

and 

Qi = T 2 {s d - s c ) = T 2 As 2 (5) 


where 


P h (E h + pP h ) - P a (E a + pP a ) cosh[/x(E„ + fiP a )/(kuT \)] 

A ‘ v i =- ^r -- + Nk B In — ' ' (6) 

T\ cosh [n(E h + fiP h ) (k B Ei) 


A Pc(E, + pp c ) - Pd(E d + pp d ) , Ar , , cosh[/i(E d + pp d )/(kn T 2 )] 

AS2= - t 2 - + Mb '" (7) 


The parameters Pj and E\ are the ferroelectric polarizations and the electric field 
strengths in the states i (= a, b, c, d) shown in Fig. 1. 



Fig. 1. The P-T diagram of a ferroelectric Stirling refrigeration-cycle. 
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From Eq. (1), one can derive 

(E + PF)/T= F(P), (8) 

where F(P) is a function only of the polarization and determined by Eq. (1). For the 
process of a given polarization P, F(P) is a constant, so that 

(E+PP)/T= const. (9) 

According to Eq. (9) and Fig. 1, one can find without difficulty that A.s'i = As 2 — 
As is a constant, which is independent of the electric field strength E and the tem¬ 
perature T. Thus, only if the change of the electric field strength with the tempera¬ 
ture is controlled properly, can the constant electric polarization process be realized. 
For two constant electric polarization processes, it has been shown [10] that 

Q\2 = Q 34, 10 

which may be derived from Eq. (3). This shows clearly that a ferroelectric Stirling 
refrigeration-cycle may possess the condition of perfect regeneration and can attain 
perfect regeneration through the use of a perfect regenerator. 

When the irreversibility of heat transfer is taken into account, the cycle is irrever¬ 
sible. The temperatures T\ and T 2 of the working substance in the high- and low- 
temperature isothermal processes are, respectively, different from the temperatures 
Th and T c of the heat sink and the cooled space, but Eqs. (4) and (5) are still true. 
When the heat transfer obeys a linear heat-transfer law [16-20], Q\ and Q 2 may be 
expressed as 


Q\ — k\(T\ - Th)t\ (11) 

and 

Qi = k 2 (T c - T 2 )t 2 , (12) 


where k\ and k 2 are the thermal conductances between the working substance and 
the heat reservoirs at temperatures 7), and T c , and t\ and t 2 are the times taken for 
the two isothermal processes at temperatures T\ and T 2 . 

For an irreversible ferroelectric Stirling refrigeration cycle, not only the isothermal 
processes but also the regenerative processes are affected by thermal resistances, so 
that the cycle cannot attain the perfect regeneration. Equally, it may be assumed 
that the regenerative losses Q, per cycle are given by 


Q, = (1 - rj) 


CpdT, 


Ti 


(13) 


where is the efficiency of the regenerator. In general, C P is closely dependent on the 
temperature. In order to simplify the calculation and analyze the main characteristics 
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of regenerative influence, C P may be substituted by the average value C P a of the 
specific heats at constant electric polarization between T\ and T 2 , so that Eq. (13) 
may be simply written as 

Qr = (1 - r])C PA (T x - T 2 ). (14) 

Because a ferroelectric Stirling refrigeration cycle is only operated between two 
heat reservoirs at temperatures 7/, and T c , the regenerative losses O, must be com¬ 
pensated from the heat sink during one constant electrical polarization process and 
released to the cooled space during the other constant electrical polarization process. 
If not, the regenerator would not work normally. 

Because the ferroelectric Stirling refrigeration-cycle is operated between two heat 
reservoirs, there naturally exists a heat leak from the heat sink at temperature 7), to 
the cooled space at temperature T c . It is often assumed [21-23] that the heat-leak 
loss is directly proportional to the temperature difference between the two heat 
reservoirs and the cycle time, i.e. 


Ql = k L {T/, — T c )t 


(15) 


where k p is the coefficient of the heat leak and t is the cycle period. 

Thus, the net amounts of heat released to the heat sink and absorbed from the 
cold reservoir per cycle are 

Qh = Qi- Q, - Ql (16) 

and 

Qc = Qi~ Qr - Ql- (17) 

Owing to the influence of thermal resistances, the regenerative processes need a 
non-negligible time compared with the time of the isothermal processes. It is often 
assumed that the time of the regenerative processes is directly proportional to the 
temperature difference of the working substance in the two isothermal processes, i.e., 

h = h = u(Ti - T 2 ), (18) 

where u is a proportionality constant which is independent of temperature. There¬ 
fore, the cycle period is given by 

t — t\ + t 2 + 2t 3 (19) 

Such a cycle model includes the influence of thermal resistances between the 
working substance and the heat reservoirs, the regenerative losses in the regenerative 
processes and the heat leak loss between the heat reservoirs. It may be used to dis¬ 
cuss the influence of several key irreversibilities on the performance of a ferroelectric 
Stirling cycle and obtain some optimal results. 
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4. General performance characteristics 

Using the above equations, one can obtain the expressions of the important per¬ 
formance parameters of the ferroelectric Stirling refrigeration cycle. For example, 
the coefficient of performance, cooling rate, and power input are, respectively, given 
by 


Qc _Qi~ Qr — Ql _ t 2 — a 
Qh Qc Q\ — Qi T\ — t 2 

Qc T 2 - b(T x - t 2 ) 

t ~ Ty , T 2 qL ' 

ki(Ti - T h ) + k 2 (T c - T 2 ) + l( 1 2) 


( 20 ) 

( 21 ) 


and 


P = 


Qh - Qc 


Ti 


+ - 


Ti 


k\{T\ - T h ) k 2 (T c - T 2 ) 


■ + b\(T\ - T 2 ) 


( 22 ) 


where a — Q L /As, b — (1 — r])C PA /As, b\ — 2m/ As, and q L — k L (T h — T c ) is the rate 
of the heat leak. 

Using a refrigerator, one always wants to obtain a cooling rate as large as possible 
for a given power input. For this purpose, we introduce the Lagrangian 


L — r + vp = 


1 + (v — b)(y — 1) 


y 


1 


- 9l, 


k\{yT 2 — Th) k 2 (T c — T 2 ) 


+ bi(y- 1) 


where v is a Lagrange multiplier and y — T\/T 2 . 
Using the extremal condition 


(23) 



(24) 


and Eq. (23), one can obtain an optimal relation as 


y(T c - T 2 ) = h(yT 2 - T h ), (25) 

where h — cjk\/k 2 . Substituting Eq. (25) into Eqs. (20)-(22) and eliminating T 2 , one 
has 
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_ T c y + l:(T h - ay) - ay 
0 DC/: r • /;7 ;i ) 

r = kT __ 

c y /(y-T h /T c ) + e(y- 1) 


(26) 

(27) 


and 


p = kT c 


_ 1 _ 

y/{y-T h /T c ) + e(y-\y 


(28) 


where A’ = 1/(1 /\fk[ + 1 /y/ki) 2 and e — b\kT c . 

Using Eqs. (26)-(28), we can generate the r* — P*, s — p* and s — r* curves, 
as shown in Figs. 2-4, respectively, where r* — r/kT/, and p*—p/kTh are the 
dimensionless cooling rate and power input. 

It is seen from Fig. 2 that the cooling rate is not, in general, a monotonic function 
of the power input, so there exists a maximum value for the cooling rate. It can be 
derived from Eq. (27) that when 


d+ e 
b + e 


T h /T c , 


(29) 



Fig. 2. The dimensionless cooling rate r* versus power input p*. Curves a (c/l = 0, b = 0, e = 0), b 
(q L = 0, b = 0, e = 0.05), c (q L = 0, b = 0.1. e = 0), d (q L = 0, h = 0.1, e = 0.05), e (q L = 0.1, b = 0, 
e = 0). f (q L = 0.1, b = 0,e = 0.05), g (q L = 0.1, b = 0.1. e = 0). and g (q L = 0.1, b = 0.1, e = 0.05) are 
plotted for T/,/T c = 2. 
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Fig. 3. The COP e versus power input p* curves. Curves a (a/7;, = 0, b = 0, e = 0), b (a/7;, 
e = 0.05), c (a/T), = 0, h = 0.1, e = 0), d (a/7 - ;, = 0, h = 0.1, e = 0.05), e (a/7* = 0.05, h = 
(a/7;, = 0.05, b = 0, e = 0.05), g (a/7;, = 0.05, b = 0.1, e = 0). and g (a/7;, = 0.05, b = 0.1. e 
plotted for 7*/7 c = 2. 



= 0, ft = 0, 
i, e = 0), f 
= 0.05) are 


Fig. 4. The COP e versus dimensionless cooling rate r* curves. The values of the parameters T;,/7 C , a/7;„ 
ft and e are the same as those used in Fig. 3. Curves a, b, c, d and e, f, g. h correspond to the cases of = 0 
and 0.1, respectively. 
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the cooling rate attains its maximum value, i.e., 


(b + e)T c 


= kT r 


(d + e)T h - (b + e)T c 


(b + e)(d + e)T c 
(d - b)[(d + e)T h - (b 4 


- Qli 


e)T c \ 


+ e 


(30) 


where d — ^(1 + b)(b + e)T c /Ti, — be. In such a case, the corresponding power-input 
and coefficient of performance are 


_ kTc _ 

(b + e)(d+e)T c 

(d - b)[(d + e)T h - (b + e)T c ] +e 


(31) 


and 

(b 4- e)T c — a(\ + h)(b + e)(d + e)/[d + e + h(b 4- e)\ 
(d+ e)T h - (b + e)T c 


(32) 


and the corresponding temperatures of the working substance in the two isothermal 
processes are 


T\ 


d 4- e 4- h(b + e) _ 
(l + h)(b + e) 


and 


T 2 


d 4- e 4- h(b + e) 
c (1 + h)(d + e) 


= T lr 


(33) 


(34) 


It is seen from Fig. 3 that the coefficient of performance is a monotonically 
decreasing function of the power input. When the power input is equal to p e , the 
coefficient of performance reduces to 0. From Eqs. (26) and (28), one obtains an 
explicit expression for p e as 


_ rj-i y £ A 

Pe c y e /i}’s - Th/T c ) + e(j’e - 1) 


(35) 


where 

y e — {(14- b)T c - a( 1 + h) - bhT h 

+ fyl + b)T e - a{ 1 4- h) - bhT h f+ 8(1 + b)bhT h T7\/{2bT c ). 


Thus, the region of p ^ p E is not allowed for a ferroelectric Stirling refrigeration 
cycle. 
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It is seen from Fig. 4, that when the regenerative losses and regenerative time are 
negligible, b — 0 and b\ — 0. The coefficient of performance is a monotonically 
decreasing function of the cooling rate. This is similar to the case of a Carnot 
refrigeration cycle, in which the time of the adiabatic processes is negligible com¬ 
pared with that of the isothermal processes. For other cases, the coefficient of per¬ 
formance is not a monotonically decreasing function of the cooling rate. There exists 
a maximum cooling rate, which is given by Eq. (30). When r < r max , there are two 
different coefficients-of-performance for a given cooling rate, where one is smaller 
than e,. and the other is larger than e r . Obviously, the region of e < s r is not optimal 
because the coefficient of performance decreases with the decrease of the cooling 
rate. For a ferroelectric Stirling refrigeration-cycle, the optimal region of the 
coefficient of performance should be determined by 

s ^ (36) 

According to Eq. (36), we may further determine the suitable ranges of the power 
input and the temperatures of the working substance in the two isothermal processes 
as 


P <Pr, 

(37) 

T1 < Ti r , 

(38) 

and 


To ^ T lr . 

(39) 

Obviously, the parameters r max , e r , p r , T [r and T 2r are very 
provide some criteria for the optimal design and operation of a 

important and may 
ferroelectric Stirling 


refrigeration cycle. 


5. Discussion 


5.1. Rapid regenerative processes 


When u = 0, this implies that the time of the two regenerative processes is negli¬ 
gible compared with that of the two isothermal processes, although the regenerative 
losses are taken into account. In this case, the r* — P*, e — p* and e — r* curves are 
shown as the curves labelled g in Figs. 2-4 and the above equations may be simplified. 
For example, Eqs. (27) (34) may be written as 


r — kT c 


\-b(y-l) 
yl(y - T h /T c ) 


~<1L, 


(40) 
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p — kT c 


y — i 

y/(y - T h /T c ) ’ 


(41) 


y = 7(1 +b)T h /(bT c ), (42) 

'max = kT c (VTTb - y/bT h /T c ) 2 - q L , (43) 

Pr = k[y/(l + b)T c /(bT h ) - 1][7), - 7WV(TT7)], (44) 

_ Ty _ a(l+/Q V(1 + \/b)T c /T h 

£r 7(1 + 1 /b)T h T c - T c 7(1 + \/b)T h T c - T c 7(1 + l/b)T c /T h + h 

- b, (45) 


T\r — Th 


7(1 + b)T c /(bT h ) + h 
1 +h 


and 


(46) 


T lr = T c 


1 + hy/bfdK} + b)T e ] 
\+h 


(47) 


5.2. Highly efficient regenerative processes 


When ;/ —»■ 1, the regenerative losses in the two constant electric polarization 
processes are negligible. In such a case, the r* — P*, e — p* and e — r* curves are 
shown as the curves labelled f in Figs. 2-4 and Eqs. (26), (27) and (29)-(34) may be 
simplified as 


T c y + h(T/, - ay) - ay 
(y-\)(T c y + hT h ) 


■ = kT r 


1 


y/(y - T h /T c ) + e(y - 1) 


- d/,. 


(49) 


y = ( 1 + d/e)T h /T c , 

kT c 

F max — /- — (jL 

l — e + eTh/T c + 2,/el\]T C 


/T . [1 + yT c /(eT h )\T h /T c - 1 

Pr = A- T, ----- ,.- , 

1 — e + eTiJT c + 27 eTh/T c 


(50) 

(51) 

(52) 
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T c - aj 1 + m + y/T c /{eT h )]/[\ + h + fiTJjeT,,)] 
WT c /(eT h ) + i]r,, - 


( 53 ) 


T\r = 7), 


1 + A + fiT c /(eT h ) 
1 +/? 


(54) 


and 


^ _ T 1 + h + 4T c /(eT„) 

2r C (1 + h)[l + jT c /(eT h )]' 


(55) 


It should be pointed out that such a ferroelectric Stirling refrigeration-cycle with 
perfect regeneration has the same performance characteristics as an endoreversible 
ferroelectric Carnot refrigeration-cycle in which the time of the adiabatic processes 
is also determined by Eq. (18). Thus, the optimal performance of an endoreversible 
ferroelectric Carnot refrigeration cycle can be directly derived from the above 
results. 


5.3. Rapid and highly efficient regenerative processes 


When u — 0 and q = 1, the r* — P *, e — p* and s — r* curves are shown as the 
curves labelled e in Figs. 2-A and Eq. (40) or (49) may be further simplified as 


y - E 1 _ 

' c y/{y - T h /T c ) 


~ dr.- 


(56) 


It is seen from curve e in Figs. 3 and 4 that the coefficient of performance is a 
monotonically decreasing function of not only the power input but also the cooling 
rate. When y —> 00 , P -* 00 , e = 0, and r — kT c — qi. It is clear that the cooling 
rate of a refrigerator is smaller than kT c . 

5.4. Zero heat-leak 


When Qi = 0, the r* — P*, e — p* and s — r* curves are shown as the curves 
labelled d in Figs. 2^4 and Eqs. (26), (27), (29), (30) and (32) may be simplified as 


s = 


1 


-b. 


(57) 


r = kT c 


1 - b{y - 1) 

y/(y ~ T h /T c ) + e(y - 1) ’ 


(58) 
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V max 


(b + e)T c _ b 

(d + e)T h - (b + e)T c 
(b + e)(d+e)T c 

(d-b)[(d+e)T h -(b + e)T c ] + e 


(59) 


and 


s r 


(b + e)T c 

(d + e)T h - (b + e)T c 


(60) 


5.5. Other expressions of the polarization 


When y6 = 0, Eq. (1) may be simplified as 


P — iV/itanh 


pE 

kft' 


(61) 


fji(E ~\~ (5P') piE 

At high temperatures, --——— << 1 and -—— << 1. Eqs. (1) and (61) return to 

k B T k b T 

the Curie-Weiss law [24-27] 


P = 


CE 

T-T 0 


(62) 


and the Curie Law [5,24,25,27] 

P = CE/T. (63) 

respectively, where To = Np 2 p/k B is the Curie temperature and C = Np 2 /k B is a 
constant characteristic of the ferroelectric materials called the Curie constant. This 
implies the fact that when the working substance obeys the Curie-Weiss law and the 
Curie Law, the optimal performance of a ferroelectric Stirling refrigeration cycle can 
be directly derived from the results in this paper. 


6. Conclusions 

We have established an irreversible cycle model of the Stirling refrigerator using a 
ferroelectric material as the working substance, based on a general expression of the 
polarization of ferroelectric materials and a linear heat-transfer law. The influences 
of irreversibilities due to thermal resistances between the working substance and the 
heat reservoirs, regenerative losses in the two regenerative processes, and heat leak 
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loss between the heat reservoirs on the performance of the ferroelectric Stirling 
refrigeration cycle have been investigated. The maximum cooling rate and other 
relevant performance parameters are calculated. Some parameter criteria, which are 
helpful for the optimal design and operation of ferroelectric refrigerators, are given. 
The results obtained are quite general. The optimal performance of not only the 
ferroelectric Stirling refrigeration cycle but also the ferroelectric Carnot refrigeration 
cycle affected by the finite-rate heat transfer may be directly derived from the results 
in this paper as long as the dielectric materials, obeying Eq. (1), (61), (62), or (63), 
are used as the working substance. 
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